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ABSTRACT 

Sgr A* at the Galactic Genter is a puzzling source. It has a mass M = (2.5 ± 
0.4) X 1O®M0 which makes it an excellent black hole candidate. Observations 
of stellar winds and other gas flows in its vicinity suggest a mass accretion 
rate M ^ few x lO“®M0yr“^. However, such an accretion rate would imply a 
luminosity > 10^° ergs“^ if the radiative efficiency is the usual 10%, whereas 
observations indicate a bolometric luminosity < 10^^ ergs“^. The spectrum of 
Sgr A* is unusual, with emission extending over many decades of wavelength. 

We present a model of Sgr A* which is based on a two-temperature optically-thin 
advection-dominated accretion flow. The model is consistent with the estimated 
M and M, and hts the observed fluxes in the cm/mm and X-ray bands as well 
as upper limits in the sub-mm and infrared bands; the £t is less good in the 
radio below 86 GHz and in y-rays above 100 MeV. The very low luminosity of 
Sgr A* is explained naturally in the model by means of advection. Most of the 
viscously dissipated energy is advected into the central mass by the accreting 
gas, and therefore the radiative efficiency is extremely low, ~ 5 x 10“®. A critical 
element of the model is the presence of an event horizon at the center which 
swallows the advected energy. The success of the model could thus be viewed 
as confirmation that Sgr A* is a black hole. 


Subject headings: accretion, accretion disks — black holes — galaxies: nuclei 
— Galaxy: center — radiation mechanisms: bremsstrahlung, inverse Gompton, 
synchrotron — radio sources: Sgr A* 


1. Introduction 


The enigmatic radio source, Sagittarius A* (Sgr A*), has for many years been a puzzle 
(Genzel & Townes 1987; Genzel, Hollenbach & Townes 1994; Mezger, Duschl & Zylka 
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1996). The source is located at the dynamical center of the Galaxy and is presumed to 
be associated with a supermassive black hole. However, observations provide conflicting 
indications both for and against the black hole hypothesis and there is currently no model 
that explains all the observations. The paradoxical clues may be summarized briefly as 
follows: 

1. Dynamical measurements indicate a dark mass of (2.5 ± 0.4) X within the 

central 0.1 pc of the Galactic Genter (Haller et al. 1996; Eckart & Genzel 1997). This is 
believed to be the mass M of the putative supermassive black hole in Sgr A*. 

2 . Observations of stellar winds and gas flows near Sgr A*, coupled with the above estimate 
of the black hole mass, suggest that Sgr A* must accrete gas from its surrounding at a 
mass accretion rate M ^ few x lO“®M 0 yr“^ (Genzel et al. 1994). 

3. This M implies an accretion luminosity L ~ O.lMc^ > 10*^° ergs“^, assuming a nominal 
radiative efficiency of 10%. However, Sgr A* is an unusually dim source with a total lu¬ 
minosity from radio to y-rays under 10^^ ergs“^. The extremely low luminosity has been 
used to argue against Sgr A* being an accreting black hole (Goldwurm et al. 1994). 

4. For the above M and M the peak emission from a standard thin accretion disk will be 
in the near infrared (cf. Frank, King & Raine 1992). In fact, this is a generic prediction of 
any model that involves an optically thick flow radiating as a blackbody. However, Menten 
et al. (1997) have obtained a strong upper limit on the 2.2 micron flux of Sgr A* that 
effectively rules out such models. 

5. Sgr A* is brightest in the radio/mm band (see §2 and Fig. 1), it is weakly detected in 
X-rays (Predehl & Trumper 1994), and it may have been detected between 100 MeV and 
2 GeV (Merck et al. 1996). These observations, combined with the Menten et al. infrared 
upper limit, imply a spectral distribution completely unlike a blackbody or even a sum of 
power-law components. The observations suggest that Sgr A* is optically thin, and that 
many different radiation processes may be in operation. 

Over the years a number of models have been proposed for Sgr A*. Some of these have 
been phenomenological approaches which aim to explain the radio and infrared spectrum 
without including any detailed dynamics (e.g. Falcke 1996, Duschl & Lesch 1994, Beckert 
& Duschl 1997). These models generally require mass accretion rates lower than those 
indicated by the observations. Other models (Melia 1992, 1994; Mastichiadis & Ozernoy 
1994) have attempted to incorporate dynamics, but with simplifying assumptions such 
as ignoring the angular momentum of the accreting gas. One model of Sgr A* that has 
attempted a self-consistent treatment of both viscous hydrodynamics and radiation processes 
is the advection-dominated accretion flow (ADAF) model of Narayan, Yi & Mahadevan 
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(1995). An early and qualitative discussion of this model was presented by Rees (1982). 
The present paper is based on the ADAF model. 

An advection-dominated accretion flow is one in which most of the energy released by 
viscous dissipation is stored in the gas and advected to the center, and only a small fraction 
of the energy is radiated (Narayan & Yi 1994, 1995a, 1995b; Chen et al. 1995; see Narayan 
1997 for a recent review). Most current work on ADAFs has been concerned with a branch 
of low M solutions (Ichimaru 1977; Rees et al. 1982; Narayan & Yi 1995b; Abramowicz et 
al. 1995) which is present for mass accretion rates below a few percent of the Eddington 
rate (Narayan 1997; Esin, McClintock & Narayan 1997). These low-M solutions make use 
of the standard a viscosity and assume a two-temperature plasma with an equipartition 
magnetic held. The two-temperature paradigm was introduced to astrophysics by Shapiro, 
Lightman & Eardley (1976; see also Phinney 1981 and Rees et al. 1982) and is based on the 
following two assumptions: (1) It is assumed that the bulk of the viscous energy is deposited 
in the ions, with only a small fraction of the energy going directly into the electrons. (2) 
Energy transfer from the ions to the electrons is assumed to occur via Coulomb collisions, 
with no signihcant nonthermal coupling being present (see Phinney 1981; Begelman & 
Chiueh 1988). In the low-M two-temperature ADAF model the ions achieve a nearly virial 
temperature, Tj ~ 10^^ K/''", where r is the radius in Schwarzschild units, while the electron 
temperature saturates at around 10^ — 10^'^ K for r < 100 (Narayan & Yi 1995b). The high 
ion temperature causes the gas to take up a quasi-spherical shape (Narayan & Yi 1995a), 
the gas is optically thin, and the flow does not suffer from any serious thermal or viscous 
instabilities (Narayan & Yi 1995b; Abramowicz et al. 1995; Kato, Abramowicz & Chen 
1996). 

The low-M branch of ADAFs has several appealing features for explaining the ob¬ 
servations of Sgr A* and other low-luminosity accretion systems. First, since most of the 
energy is advected with the gas and lost into the black hole, the model naturally explains 
the low luminosity of Sgr A*. Second, being optically thin, the spectrum is quite different 
from blackbody. Indeed, the high electron and ion temperatures allow a variety of radiation 
processes to operate: synchrotron, bremsstrahlung and inverse Compton from the electrons 
(Narayan & Yi 1995b; Mahadevan 1997), and y-ray emission from the ions via pion pro¬ 
duction (Mahadevan, Narayan & Krolik 1997). Therefore, the gas is likely to radiate over 
a wide range of wavelengths. 

The Narayan et al. (1995) ADAF model of Sgr A* gave a reasonable £t to the obser¬ 
vations available at that time and provided an explanation for the low luminosity of the 
source. However, the model required a black hole mass of 7 x IO^Mq which differs from 
the most recent dynamical mass estimate of 2.5 x IO^Mq (Eckart & Genzel 1997). This 
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problem prompted us to take another look at the model. 

In addition, there have been important developments on the observational and theoret¬ 
ical front which further motivate the present study. The recent Menten et ah (1997) limit 
on the infrared flux of Sgr A* is well below previous “detections” (Rosa et al. 1991; Eckart 
et al. 1992) which were used by Narayan et al. (1995) to £t the ADAF model. Similarly, 
ASCA observations by Koyama et al. (1996) have shown that the X-ray emission from the 
Galactic Center region is dominated by an X-ray burster which lies within 1.3 arcminutes 
of Sgr A* . This implies that most previous measurements of the X-ray flux of Sgr A* 
(Pavlinskii, Grebenev, & Sunyaev 1992, 1994; Skinner et al. 1987) are suspect since they 
were obtained with inadequate spatial resolution. Only the ROSAT PSPC detection by 
Predehl & Trumper (1994) survives as a reliable X-ray detection, although the correspond¬ 
ing luminosity of the source in the ROSAT band is uncertain due to uncertainties in the 
absorbing hydrogen column (cf. §2). 

The modeling techniques have also advanced significantly during the intervening two 
years. The Narayan et al. (1995) model was a local one where the dynamics were calculated 
using a local self-similar solution (Narayan & Yi 1994) and the radiation was computed with 
a local approximation for Compton scattering (Dermer, Liang & Canfield 1991; Narayan & 
Yi 1995b). Consistent global dynamical solutions have since been calculated, initially with 
a pseudo-Newtonian potential (Narayan, Kato & Honma 1997; Chen, Abramowicz & Lasota 
1997) and more recently for the full relativistic Kerr problem (Abramowicz et al. 1996; Peitz 
& Appl 1997; Gammie & Popham 1997). A global scheme for Comptonization has also been 
developed (Narayan, Barret & McClintock 1997), based on techniques due to Poutanen & 
Svensson (1996). Finally, Nakamura et al. (1997) have shown that in addition to energy 
advection by ions, which had been considered in previous work, advection by electrons 
can be important in some circumstances. This effect is now included in the calculations 
(Appendix A). The present models are thus superior to those used by Narayan et al. (1995). 

In this paper we present improved ADAF models of Sgr A*. We make use of the 
spectral data described in §2 and employ the modeling techniques outlined in §3. We 
present detailed results in §4 and conclude with a discussion in §5. 


2. Spectral Data 

The radio to near infrared (NIR) spectrum of Sgr A* has been of constant interest 
since the discovery of the source by Balick & Brown (1974), and observations have been 
carried out from 400 MHz (Davies et al. 1976) to ~ 10^"^ Hz (Menten et al. 1997). Table 
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1 is a compilation of the radio to NIR observations. The data are also plotted in Fig. 1. 
The two general features in the table and hgures are 1) we have given only the maximum 
and minimum fluxes for frequencies at which variability has been observed, and 2) open 
circles correspond to low resolution measurements, which we treat as upper limits, while 
hlled circles correspond to the highest resolution measurements in each band (subject to 
the additional constraint that we require the resolution to be < 1 arcmin). The latter 
measurements, which are identihed with a in Table 1, are most likely to probe the actual 
accretion flow in Sgr A*, while the former could be contaminated by emission from other 
components in the source. This is discussed in detail below. 

In determining the spectrum of the ADAF from Sgr A*, special attention has to be 
paid to the angular resolution of the observations, which ideally must be comparable to 
the size of the ADAF. The angular size of the ADAF in Sgr A* depends on the wavelength 
band under consideration. The low frequency (< 1 GHz) radio emission is from radii on 
the order of ~ — 10"^(where Rg is the Schwarzschild radius), which corresponds 

to an angular size at the Galactic Genter of ^ 0.06" (~ 500 AU), where we have taken 
M = 2.5 X IO^Mq. With increasing frequency, the size becomes smaller, and the sub-mm 
and infrared emission are from no more than a few tens of Rg. The soft X-rays are again 
from a large volume (up to about 10"^/?^), but harder X-rays are from progressively smaller 
radii. Gamma-ray emission from pion decay is again limited to a few tens of Rg. 

The radio spectrum of Sgr A* appears to consist of two components, with a break 
at around the z/ ~ 86 GHz point in Fig. 1. The spectral dependence is Ljj ~ for 
V < 40GHz, and steepens to Ljj ~ ® above 40 GHz. The steep increase has been noted 

by Zylka et ah (1992, 1995), Serabyn et ah (1992, 1997) and Rogers et al. (1995), but the 
origin of the two components is unclear. Perhaps coincidentally, for frequencies ^ 50 GHz, 
scattering by turbulent density fluctuations in the ISM leads to source broadening. The 
apparent size of the source, which is proportional to the square of the observed wavelength 
(e.g. Davies et al. 1976; Backer 1982), is larger than the intrinsic size at these frequencies. 
The low frequency data could thus in principle be contaminated by structures in the source 
(e.g. jets) which are larger than the relevant radiating region of the ADAF but smaller 
than the scatter-broadened size. 

The ADAF model we present in the following sections hts the data above 86 GHz 
well but underpredicts the flux at lower frequencies. Motivated by this, and in view of the 
scattering argument mentioned above, we plot the low frequency data as open circles in 
Fig. 1. We treat these data as upper bounds rather than as hrm detections of the central 
accretion flow in Sgr A*. Thus, even though the Marcaide et al. (1992) and Alberdi et al. 
(1993) measurements at 22 GHz had exceptionally good resolution, we still plot these points 



- 6 - 

as open circles because they are affected by scattering. At 22 GHz, the scattered source 
size is ~ 15AU, whereas the relevant region of the ADAF which produces the radiation has 
a size of ~ lAU. Any contamination to the radio flux must be present on a scale ^ 15 AU, 
which is rather a stringent requirement. We must therefore keep open the possibility that 
the low frequency radiation (i.e. below the apparent break at ~ 86 GHz) also arises from 
the accretion flow. The models presented in this paper are then unable to account for the 
observed level of this low frequency emission. 

For frequencies ^ 50 GHz, the scattering size falls below the intrinsic size. High reso¬ 
lution observations can therefore determine directly the emission from the central engine. 
However, only a few high resolution 1 mas) observations have been made at these wave¬ 
lengths using VLBI or the VLBA. These radio observations are represented as filled circles 
in Fig. 1. We take the point at 86 GHz to be the best determined radio flux of Sgr A* for 
the purposes of modeling the source. Apart from the VLBI/VLBA points, all other radio 
points ^ 86 GHz have been represented by open circles since these observations were done 
with poor resolution 4"). 

The FIR to NIR observations are all upper limits, and the low angular resolution points 
are represented by open circles. The upper limits by Gezari et ah (1994) and Zylka et al. 
(1992) are the best resolution observations in the FIR band, and are therefore shown as 
filled circles. Stolovy et al. (1996) detected a source in the mid-infrared with an extinction- 
corrected flux of 100±40 mJy. However, their resolution was 0.5" and their flux may possibly 
be contaminated by the compact cluster of sources observed by Eckart et al. (1995). We 
plot their data point as an upper limit with a hlled circle. The recent upper limit of Menten 
et al. (1997) in the NIR at 2.2/rm is a very high resolution measurement since the authors 
were able to determine accurately the relative position of Sgr A* in their speckle images 
by comparing with a radio map of nearby stellar SiO masers. The effective resolution of 
the speckle image is 0.15", and therefore we plot this limit again as a hlled circle. 

X-ray and hard X-ray as well as gamma-ray observations of Sgr A* have long been 
limited by the relatively poor angular resolution available at these energies and the very 
crowded held of the Galactic Genter region. Even at hard X-ray energies (e.g. 20-100 keV) 
where the density of sources is relatively low, the Sgr A* source region contains more sources 
(~10) in the central 5° of the Galaxy than any other region (Goldwurm et al. 1994). At 
soft X-ray (e.g. 0.5-4 keV) and low-medium X-ray energies (e.g. 2-10 keV), this same 
region is even more crowded, with complex dihuse emission and at least 10 sources within 
the central 1° of Sgr A* (cf. Watson et al 1991 and Koyama et al 1996). Thus the hux 
and spectral distribution of the Galactic Genter in the ~1-100 keV band, where the ADAF 
bremsstrahlung component peaks (cf. Fig. 1), is particularly critical and angular resolution 
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is paramount. At the soft X-ray end of this band, the relatively uncertain absorbing column 
density, N^, towards the Sgr A* source is especially important since the value chosen (which 
can in principle be determined by future high resolution observations with AXAF) greatly 
affects the derived source flux and luminosity. 

We have included in the X-ray/hard X-ray region of the spectrum plotted in Figure 1 
only 4 points or upper limits which all are the highest angular resolution available in their 
respective energy bands. These data are listed in Table 2. The one detection of a point-like 
source with position fully consistent with Sgr A* is the ROSAT detection (PSPC) with ~ 20" 
resolution in the 0.8-2.5 keV band reported by Predehl and Trumper (1994). This detection 
is plotted for an assumed = 6 x 10^^ cm“^, which is the usual best estimate (cf. Watson 
et al 1981) for the interstellar column density and which corresponds to the usually quoted 
visual extinction of Ay ~ 25-30 mag. For this Nh, the total integrated luminosity in the 
0.8-2.5 keV band is Lx = 1.5 x 10^^ erg s“^. [| The used here differs from the much 
higher value (1.5 x 10^^ cm“^) assumed by Predehl and Trumper. These authors chose a 
higher Nh in order to make the soft X-ray luminosity more compatible with the variable 
compact source within ~ F of Sgr A* reported by Skinner et al (1987) and Pavlinskii, 
Grebenev and Sunyaev (1992) from 2-20 keV coded aperture imaging observations with 
modest (few arcmin) resolution. However, recent ASCA observations of the Galactic center 
in the 2-10 keV band, with ~ 1' angular resolution but much higher spectral resolution, 
have been reported by Koyama et al. (1996) and Maeda et al (1996). They show the Sgr 
A* region to be complex, with diffuse emission over a ~2 x3' region and a point source 
at 1.3' from Sgr A* which is an X-ray burster and therefore a neutron star in a (8.4h) 
binary system. The burster has very likely dominated the coded aperture imaging flux 
measurements. Koyama et al. therefore quote an upper limit of 10^® erg/s as the 2-10 keV 
luminosity for the entire Sgr A* complex, with a possible actual value of ~10^^ erg/s for 
the point source alone. We have plotted this ASGA upper limit in Fig. 1. 

The ROSAT PSPG detection, although only with angular resolution of ~ 20", is a likely 
detection of Sgr A* with minimal contamination from the surrounding diffuse source since it 
appears in the ROSAT image to be consistent with a point source. However, higher spatial 
resolution (e.g. AXAF) X-ray imaging, and/or evidence for time variability, are needed to 
conhrm this identideation. In this sense, even the ROSAT detection might be regarded as 
an upper limit. Nevertheless, this data point is plotted as a solid symbol as it corresponds 
to the highest resolution X-ray observation. The luminosity plotted in Fig. 1 and listed in 
Table 2 corresponds to the likely Nh value of 6 x 10^^ cm“^ and was obtained by considering 


^The Lx and calculations were done for the ROSAT data using the PIMMS program, supplied 
by the HEASARC at GSFC. 



two extreme values of a power law spectral index: photon index 1.0 and 2.0. These indices 
bracket the corresponding index (~1.4) for a bremsstrahlung spectrum as predicted by our 
model. The detected flux in the soft (0.8-2.5 keV) band is, however, still highly sensitive 
to the assumed Nh (as also pointed out by Predehl and Trumper): for Nh = 8 x 10^^ the 
flux increases by a factor of 2.8 over that shown; whereas for Nh = 5 x 10^^, it decreases 
by a factor of 0.6. Thus the hts to ADAF models described in the next section, which are 
otherwise greatly constrained by this ROSAT detection, should be regarded as uncertain 
in normalization by a factor perhaps as large as 2 if uncertainties in the interstellar Nh are 
allowed for completely and perhaps as large as 3 if internal self-absorption in the vicinity 
of the AGN is allowed for (as suggested by Predehl and Trumper). The box around the 
ROSAT point in Fig. 1 includes both the uncertainty due to spectral index as well as 
an additional factor of 2 uncertainty in normalizaton. The uncertain value of Nh can, in 
principle, be directly measured or at least greatly constrained by observations of Sgr A* 
with the ACTS instrument on AXAF. This could observe the interstellar absorption edge 
of Oxygen (at 0.8 keV) at the high angular resolution (1") needed to isolate completely 
the surrounding diffuse emission. 

At still higher energies, 20-100 keV, we plot the two upper limits derived from the 
deep SIGMA observations of the Galactic Genter complex reported by Goldwurm et al. 
(1994). These observations are not able to resolve the burster, and are likely affected by it 
(or, rather, constrain the hard X-ray flux from the burster). More sensitive observations 
in the hard X-ray band at energies above 60-100 keV, where the hard X-ray spectral 
component of neutron stars in bursters usually is not detectable (e.g. Barret & Grindlay 
1995) could provide a more stringent test of the ADAF spectrum as well as the predicted 
spectral turnover. 

The y-ray spectrum shown in Fig. 1 corresponds to emission from the EGRET source 
2EG 1746-2852 which is coincident with the Galactic Genter (Merck et al. 1996). The 
source is described as point-like to within the resolution of the instrument (~ 1°), and is 
a signihcant excess (~10 (t) above the local diffuse emission. The source has a very hard 
spectrum with a photon index of 1.7 which differs signihcantly from other unidentihed 
Galactic EGRET sources. The source might therefore represent emission from Sgr A*, but 
it could equally well be unresolved emission from a giant molecular cloud in the vicinity of 
the Galactic Genter. If the source is Sgr A*, the flux and spectrum may be explained as 
y-ray emission via pion production and decay in the ADAF (Mahadevan et al. 1997). If 
the source is instead associated with a dense cloud of molecular hydrogen situated at the 
Galactic Genter, it will again show a pion decay spectrum but with a spectral index (at 
energies above 100 MeV) given by the local cosmic ray proton spectrum (i.e. a power law 
with energy index ~ —1.6). Due to the low resolution of EGRET (1° beam) and uncertainty 



in the nature of the source, we plot the 7 -ray data in Fig. 1 as open circles. The EGRET 
observations (as reported by Mayer-Hasselwander, 4th Compton Symp.) indicate possible 
variability, which would of course suggest that the source is indeed Sgr A*, but the same 
EGRET data also indicate that the source may be slightly resolved, which would argue 
that the radiation is signihcantly contaminated by surrounding diffuse emission. Future 
observations with higher sensitivity and better resolution are needed to resolve this issue. 


3. Modeling Techniques 

We consider a black hole of mass M accreting through an ADAF at a rate mM-Edd, 
where MEdd = lOTEdd/c^ = 1-39 x 10^®(M/Mq) gs“^ = 2.21 x 10“®(M/Mq) An 

ADAF is in some sense dynamically intermediate between a thin disk and a spherical 
accretion flow: it is hot and quasi-spherical {H/R ~ 1, Narayan & Yi 1995a), with rapid 
radial inflow, yet centrifugal support and viscous angular momentum transport still play 
a signihcant role. In describing the ADAF we shall generally refer to the scaled radius 
r = R/Rs, where Rg = 2GM/c^ is the Schwarzschild radius. We assume that the ADAF 
extends from the black hole horizon, r = 1, to an outer radius rout- We take the mean 
angular momentum vector of the ADAF to be inclined at an angle i to the line of sight. 

To hnd the dynamical structure of the ADAF, we use the fully relativistic, self- 
consistent, steady-state global solutions in the Kerr metric developed by Popham & Gammie 
(1997). Their model uses a nearly standard viscosity prescription which is parameterized 
by the efficiency of angular momentum transport, a (assumed constant, independent of 
radius). The viscosity prescription has been modihed, however, so as to preserve causality 
and to include relativistic effects. We also use a quasi-spherical prescription for the vertical 
structure based on Abramowicz, Lanza, and Percival (1997). Gomparison of relativistic and 
nonrelativistic global models, and the self-similar solution of Narayan & Yi (1994) shows 
that they are rather similar outside about r = 5. Inside this radius the density in our 
solution is very close to self-similar, but the temperature is slightly lower. 

The dynamical model is uniquely specihed by four structure parameters: a, the viscos¬ 
ity parameter; 7, the adiabatic index of the fluid, which is assumed to be a mixture of gas 
and magnetic helds; /, the advection parameter, which gives the ratio of advected energy 
to the viscous heat input; and a, the rotation parameter of the black hole. Of these four pa¬ 
rameters, 7 is determined in terms of the magnetic field equipartition parameter (3 (defined 
below) according to (Esin 1997): 7 = (8 —3/3 )/(6 —3/3). The viscosity parameter a is again 
determined in terms of /3 following Hawley, Gammie & Balbus (1996): a = 0.6(1 —/3). The 
parameter / is determined self-consistently by iterating between the dynamical solution 
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and the radiation solution and seeking consistency between the two (see Esin et ah 1997 
for details). The models considered here are extremely advection-dominated, so that / is 
very close to unity in all cases. (For instance, for the baseline model shown in Fig. 1, we 
obtain / = 0.9994.) Finally, we consider only the case a = 0, i.e. a Schwarzschild black 
hole. 


The thermodynamic state of the flow is described by the ion temperature Tj, the 
electron temperature T^, and the magnetic pressure Pmag = (1 ~ P)'Ptot- The total pressure 
is hxed by the dynamical model. We assume that /? is a constant, independent of radius, and 
so a and 7 are also independent of radius. Thus, all the structural parameters are constant 
in our model (cf. §3). Given f3 we know the run of gas pressure with radius. Since the 
electrons can cool efficiently and coupling between the ions and electrons is assumed to be 
weak (provided solely by Coulomb collisions), we have T* 3> Tg. The run of gas pressure 
then gives Tj to a good approximation. 

The heart of the problem is now calculating the electron temperature, which determines 
the spectral properties of the flow. A detailed description of our procedure is given in Esin et 
ah (1997). Basically, the electron temperature is determined by solving an energy balance 
equation: 

ge,adv _ gie ^ ^ gvis _ grad _ g ^ 

Here is the rate at which energy in the electrons is advected inward by the flow. This 

term is discussed in greater detail in Appendix A. Our earlier work (with the exception of 
Esin et al. 1997) ignored the advective term, but recently Nakamura et ah (1997) (see also 
Mahadevan & Quataert 1997) have shown that it can be important, particularly at low mass 
accretion rates. At low m the electrons are nearly adiabatic and the compressive heating 
as the gas flows in becomes large. The other terms in equation (1) are as follows: is the 

rate at which the electrons are heated by Coulomb collisions, Q™ is the total rate of viscous 
dissipation, and is the radiative cooling. The parameter 5 describes the fraction of 
the viscous heating that goes into the electrons. In general we set 5 = 10“^ ~ me/mp. 

Because of Compton scattering, the radiative cooling term couples the flow at different 
radii. It is calculated using the iterative scattering method described in detail in Narayan, 
Barret & McClintock (1997). The ADAF is represented by a logarithmically spaced grid of 
nested shells. Within each shell all flow variables are assumed uniform for n/2 — 6 h < 0 < 
n/2 + Oh, where 9h is an effective angular scale height, calculated as in Narayan, Barret 
& McClintock (1997) Appendix A. We then guess a value of the electron temperature 
and compute the synchrotron and bremsstrahlung emission from each shell. Then, using 
the probability matrix elements Pjj, which give the probability that a photon emitted in 
shell i of the ADAF is scattered by an electron in shell j, the rate of cooling through 
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Compton scattering is calculated (see Narayan, Barret & McClintock 1997 for details). 
This procedure is iterated until convergence is achieved. The iteration method ensures that 
multiple scattering within the ADAF is properly taken into account. It also ensures that the 
advection parameter / used in the dynamical solution is consistent with the actual energy 
advection, as determined by the radiation solution. Doppler shifts and ray deflections are 
ignored in the present calculations. 

Once the cooling term is calculated, we must get the photons out to the observer. 
We use essentially Newtonian photon transport, although we have taken a step toward 
relativistic photon transport by including gravitational redshift (this is also included in the 
Compton scattering calculation). We do not include Doppler shifts, which to lowest order 
in vjc broaden the spectrum. Since all features in our spectra have dvjv of order unity, 
this should not greatly alter the gross properties of the spectrum, except where it is very 
steep. At next higher order in v/c gravitational redshift, higher order Doppler, and various 
geometric effects all enter. In our model we have included a correction for the volume of 
the emitting region (1/^1 — 1/rx the Euclidean value) and for gravitational redshift. The 

sum of these effects is that Ly(z/) observed at large radius is now — 1/r), that is, 

Li, is shifted redward by the gravitational redshift factor. Jaroszynski & Kurpiewski (1997) 
have recently considered the full effects of photon transport near a Kerr black hole on the 
spectrum of an ADAF. 

The above discussion is concerned with radiation from electrons. We also compute 
pion production by the hot protons in the ADAF and the resulting y-ray emission through 
pion decay. The procedure we follow for this calculation is described in Mahadevan et al. 
(1997) and is based on earlier work by Dermer (1986). The present calculations differ in 
three respects from those in Mahadevan et ah: 1) the density prohle of the protons is now 
given by the relativistic global solution of Popham & Gammie (1997) instead of the pseudo- 
Newtonian solution of Narayan, Kato & Honma (1997), 2) we allow for the non-spherical 
geometry of the flow via the angle Oh mentioned above, and 3) we include the effect of 
gravitational redshift. 

The pion production rate depends sensitively on the energy distribution of the protons. 
Mahadevan et al. (1997) considered two extreme distributions: a thermal distribution and 
a power-law distribution. The latter has signihcantly more pion production and maximizes 
the efficiency of y-ray emission. We employ the power-law model with an energy index 
p = 2.3. The choice of model is dictated by the EGRET spectrum of the Galactic Genter, 
which is consistent with a power-law distribution of proton energies and is inconsistent 
with a thermal distribution (Mahadevan et al. 1997). Note that Mahadevan & Quataert 
(1997) have shown that the protons in an ADAF do not have time to thermalize; therefore. 



it is permissible to assume a nonthermal energy distribution for the protons. 


4. Results 

In this section we combine the modeling techniques outlined in §3 with the observa¬ 
tional constraints described in §2 to come up with ADAF models of Sgr A*. Our standard 
parameters are as follows. We assume that the ADAF extends from r = 1 to rout = 10^. 
For the black hole mass, we take the estimate of Eckart & Genzel (1997): M = 2.5 x lO^M©. 
We assume exact equipartition between gas and magnetic pressure in the accreting gas: 
(3 = 0.5. Following Hawley, Gammie & Balbus (1996) we take the viscosity parameter to 
be a ~ 0.6(1 — {3) = 0.3 (the coefficient lies in the range 0.5 — 0.6). We assume that elec¬ 
trons receive only a small fraction 5 ~ me/mp of the viscous dissipation as direct heating: 
6 = 0.001. We assume a Schwarzschild black hole: a = 0. Since we do not know the ori¬ 
entation of the angular momentum vector of the accreting gas, we set the inclination of 
the system to a generic value: i = 60°. The models considered here are extremely optically 
thin and so the results are virtually independent of i. 

The only parameter that needs discussion is rout- In a viscous rotating accretion flow, 
the outer edge of the flow is determined by the radius at which the outward viscous transport 
of angular momentum is balanced by external torques. In the case of Sgr A* it is not clear 
exactly where this balance is achieved, but it is likely to be roughly at the radius where 
the mass supply originates. Our choice of rout = 10^, corresponding to i?out = 0.024 pc, 
may be an underestimate of the true outer radius. However, this is not a serious concern 
for the modeling since for such large outer radii the results are insensitive to the actual 
value chosen. This is because the energetically important region of the flow is closer to the 
black hole. 

As described above, all the parameters of the model are fixed at standard values in 
our model. The only parameter that we consider fully adjustable is the mass accretion rate 
rh. In each model described below, we have adjusted rh so as to £t the ROSAT X-ray 
flux (but see the discussion in §2 for uncertainties in the X-ray flux due to the hydrogen 
column). Although we treat m as a free parameter, there are in fact some constraints. The 
infrared source IRS 16 appears to be the primary supplier of gas to Sgr A*. Assuming 
that Sgr A* captures a fraction of the wind of IRS 16 by Bondi accretion, Melia (1992) 
estimated M ^ 2 x 10“^ assuming a wind velocity of 600 kms“^, while Genzel et 

ah (1994) estimated a mass accretion rate of M ^ 6 x 10“® using a wind velocity 

of 1000 kms“^ . We take these two estimates to be low and high extremes. Converting to 
Eddington units, m is thus likely to be in the range 10“^ < m < 3 x 10“^. 
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The solid line in Fig. 1 shows our baseline model where M, a, (3, and 6 have been 
set to their standard values (2.5 x IO^Mq, 0.3, 0.5, 0.001 respectively) and rh has been 
adjusted to £t the X-ray flux. This model has a mass accretion rate of m = 1.3 x 10“^, which 
lies within the acceptable range of m discussed above. The computed spectrum has four 
well-defined peaks. From the left, these correspond to self-absorbed thermal synchrotron 
emission, singly Compton-scattered synchrotron radiation, bremsstrahlung emission, and 
7 -rays from neutral pion decay. The model spectrum passes through the VLBI radio flux 
measurement at 86 GHz, which we have identified to be a reliable high frequency radio 
observation (cf. §2), and turns over sharply in the sub-millimeter band exactly as required 
by the sub-millimeter and infrared upper limits. The model also satisfies the ROSAT X-ray 
detection and the other X-ray upper limits. Considering that only one parameter, m, has 
been adjusted, we consider the agreement with the data satisfactory. 

Note that the infrared and X-ray fluxes have come down significantly compared to the 
data shown in Narayan et al. (1995). The present model is compatible with both the new 
measurements. Interestingly, we cannot find any model that fits either the new infrared 
upper limit with the old X-ray data, or the new X-ray flux with the old infrared data. 

The model does not agree with the low frequency radio measurements at frequencies 
below 86 GHz. We discuss this discrepancy at the end of this section. 

The model also has a discrepancy in the y-ray band; the predicted flux is lower than 
the observations by approximately an order of magnitude. The pion decay model requires 
an accretion rate of rh = 4.5 x 10“^ to fit the observed flux whereas the fit to the rest of 
the spectrum gives rh = 1.3 x 10“^. There is thus a discrepancy of a factor ~ 3 —4 between 
the two values of rh. We have been unable to come up with a reasonable resolution of this 
discrepancy. Considering that we have used an extreme model which maximizes the pion 
production rate (cf. §2), we are compelled to suggest that perhaps the EGRET source does 
not correspond to Sgr A*, but rather to a concentration of molecular gas at the Galactic 
Center. The resulting emission is then expected to be much more diffuse, but it would still 
appear unresolved to the 1° beam of EGRET. 

Two features of the ADAF model in Fig. 1 should be highlighted. First, the model 
fits the data using a reasonable mass for the black hole, M = 2.5 x IO^Mq. This is an 
improvement over the model described in Narayan et al. (1995) where the data could be 
fitted only with a mass M = 7 x IO^Mq. The primary reason for the improvement is the 
inclusion of electron energy advection, or compressive heating of the electrons (Nakamura et 
al. 1997), as we explain below. Second, the model is extremely advection-dominated. The 
bolometric luminosity Lboi integrated over all frequencies is only 2.1 x 10^® ergs“^, which 
corresponds to a radiative efficiency, e = Lboi/Afc^ = 5 x 10“®. It is this extraordinarily 
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low radiative efficiency that allows the model to fit the observations with such a large mass 
accretion rate. 

In contrast, a standard thin accretion disk model of Sgr A* is ruled out quite compre¬ 
hensively by the radio-IR-X-ray data. The dotted line in Fig. 1 shows a thin disk model 
with 771 = 10“"^, the lowest mass accretion rate that we consider reasonable. The spectrum 
was calculated assuming that the emission is blackbody at each radius. The model predicts 
an infrared flux which is many orders of magnitude above the measured upper limits. The 
short dashed line in Fig. 1 shows another thin disk model where tti has been reduced to 
10“®. This model does satisfy the IR upper limits but does not £t any of the observations, 
either in the radio or X-ray bands. Further, the mass accretion rate is unreasonably small. 

A number of improvements have taken place in the modeling techniques since the pub¬ 
lication of our first model of Sgr A* (Narayan et ah 1995). It is interesting to investigate 
what effect each improvement has had on the calculated spectrum. Figure 2 shows a se¬ 
quence of models in which we start with the simplest version of the model and progressively 
add features one by one. Figures 3-5 explore the effects of varying model parameters and 
are described below. 

The dotted line in Fig. 2a corresponds to the most primitive version of the ADAF 
model, in which the flow is assumed to have a self-similar form (Narayan & Yi 1994) 
and neither compressive heating of electrons nor gravitational redshift is included. The 
optimized accretion rate is m = 6 x 10“®. The general shape of the spectrum is similar 
to that of our standard model (the solid line), but this model differs in three ways. First, 
the emission is stronger in nearly all bands compared to our standard model; in fact, the 
model is inconsistent with the infrared limit. Second, the mass accretion rate is lower than 
in our standard model by a factor of 2. This, combined with the higher luminosity, means 
that the model is not as advection-dominated as our standard model; we obtain a radiative 
efficiency of e = 2.4 x 10“"^. Third, the calculated synchrotron emission is well above the 86 
GHz VLBI point. The discrepancy in the radio flux is more than an order of magnitude 
and is serious. Within the framework of the self-similar flow assumption, the only way to 
eliminate this problem is by changing one or more of the model parameters. This was, in 
fact, the primary reason why the Narayan et ah (1995) model required a black hole mass 
of 7 X IO^Mq instead of the current best estimate of 2.5 x IO^Mq (see Fig. 3 below which 
shows how the spectrum is modified when the black hole mass is reduced). 

The short-dashed line in Fig. 2a shows the effect of replacing the self-similar flow by a 
global flow based on a pseudo-Newtonian potential (Narayan, Kato & Houma 1997; Chen, 
Abramowicz & Lasota 1997). The changes are minimal. 
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The long-dashed line in Fig. 2a next shows the effect of including compressive heating 
(Appendix A and Nakamura et al. 1997). This model, which has rh = 1.1 x 10“^, has 
a very different spectrum compared to the previous two models. The overall emission is 
significantly reduced, notably in the synchrotron peak, and the flow is substantially more 
advection-dominated: e = 2.1 x 10“®. This model fits the 86 GHz data point very well 
without requiring any adjustment to the black hole mass. We thus confirm the result of 
Nakamura et al. (1997) that compressive heating is an important effect. Figure 2b shows 
the electron temperature profiles of the various models. The model with compressive heating 
has a significantly different temperature structure than the models without. At very low m, 
as in our model of Sgr A*, the dominant terms in the electron energy equation are the two 
pieces of written in the Appendix, viz. the terms proportional to dTe/dR and dp/dR. 

These two dominant terms balance each other, while the rest of the terms in equation (1) 
are small. In other words, the electrons are essentially adiabatic. The adiabatic condition 
gives a lower electron temperature than in the previous two models and this accounts for 
the change in the shape of the synchrotron contribution to the spectrum. 

The dash-dot line in Fig. 2a shows next the effect of including a fully relativistic global 
solution (rh = 1.3 x 10“^), taken from Popham & Gammie (1997), instead of the pseudo- 
Newtonian global solution employed in the previous two models. Most of the features are 
similar, but the overall emission is increased because the relativistic solution has a lower 
radial velocity close to the black hole compared to the pseudo-Newtonian solution. (This 
is because v < c at all radii in the relativistic model whereas the pseudo-Newtonian model 
gives V > c close to the black hole.) Gonsequently, the density is higher in the relativistic 
model and this leads to increased emission. 

Finally, the solid line in Fig. 2a shows our standard model, which includes gravitational 
redshift. As expected, this model has a lower luminosity than the previous model, but is 
otherwise quite similar. 

Figures 3-5 show the effect of varying the other parameters in the model. Figure 3a 
shows what happens when we vary the mass. In addition to the baseline model, we present 
two other models with M = 1/3 and 3 times M of the baseline model. As Mahadevan 
(1997) has shown, with increasing M, the magnetic held strength decreases as and 

this causes the synchrotron peak to move to lower frequencies oc Of the three models 

shown in Fig. 3a, the one with M = 2.5 x IO^Mq gives the best fit to the VLBI radio data. 
Assuming our choices of a and (3 are correct, this provides additional support for the black 
hole mass measured by Haller et al. (1996) and Eckart & Genzel (1997). 

Figure 3b shows the effect of varying the mass accretion rate. In addition to the 
baseline model, four other models are shown, with m = 1/2, 1/v^, and 2 times m 



of the baseline model. The radiative efficiency of ADAFs varies rapidly with m: e oc m 
(Narayan & Yi 1995b). Therefore, the luminosity varies as tv?. This can be seen in both 
the synchrotron and bremsstrahlung peaks. The Compton peak in the infrared shows an 
even stronger dependence on m. This is because the Comptonized flux is proportional to 
the product of the synchrotron emission (which is oc w?) and the optical depth (which is 
oc rh), so that the amplitude of this peak varies approximately as rh^. The three intermediate 
models in Fig. 3b are consistent with all the data, but the two extreme models lie outside 
the X-ray error box and one of them also violates the sub-millimeter and infrared limits. 
Since the X-ray luminosity of Sgr A* is uncertain. Fig. 3b indicates how the model m will 
need to be modihed if a future determination of Nh leads to a signihcant revision of the 
X-ray flux. 

Figure 4a shows the effect of varying a. For each a, we have adjusted the accretion rate 
to £t the X-ray flux. The spectrum is not very sensitive to a (except in the optical/IR), 
though the fitted mass accretion rates show modest variations: m = 6.3 x 10“^, 1.0 x 
10“^, 1.3 X 10“^, 1.5 X 10“^ for a = 0.1, 0.2, 0.3, 0.4 respectively. 

Figure 4b shows the effect of varying j3. Again, for each (3, we have adjusted the 
accretion rate to fit the X-ray flux. In this case, we see quite substantial changes in the 
predicted spectrum. An increase in (3 leads to a decrease in the magnetic held strength 
and one may be tempted to think that this would cause the synchrotron peak to reduce in 
amplitude. In fact, the opposite behavior is seen. The reason can be traced to compressive 
heating. Since the electrons are effectively adiabatic, their temperature profile is determined 
by their adiabatic index 7 e. We show in Appendix A that because the gas is a mixture of 
particles and magnetic held, the effective 7 e depends on the parameter (3. As f3 goes up, 
the magnetic pressure goes down and 7 e increases. A larger 7 e leads to hotter electrons, 
see Fig. 4c, which causes more synchrotron emission. Unfortunately, this means that the 
results are sensitive to the details of how we model energy advection in the electrons. In 
this context we note that our equation of state for the electrons differs somewhat from the 
one used by Nakamura et ah (1997). We use a relativistic equation of state which causes 
7 e to change as the electron temperature approaches and crosses the value Te ~ rrieC^/k. 
In addition, we modify the adiabatic index of the electrons to allow for the equipartition 
magnetic held which is assumed to be present in the gas and is coupled to the electrons. 
Nakamura et ah (1997) treat the particles as a separate component and ignore the held in 
their equation of state. 

Figure 4d shows the effect of varying 6. For all S < 0.01 the spectrum is generally 
unaffected. This result is different from that found in Narayan, Barret & McClintock (1997). 
Once again, the reason is the inclusion of compressive heating. Because the electrons are 
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now essentially adiabatic, their temperature profile is not affected by modest changes in the 
heating or cooling. Only when S is large, e.g. S = 0.03162 (the dot-dashed line in Fig. 4d), 
do the electrons experience signihcant additional heating and only then does the spectrum 
show a noticeable change. 

We turn now to the low frequency radio data, where the model deviates substantially 
from the measurements. As discussed in §2, this region of the spectrum may possibly arise 
from a separate component such as a jet which is outside the scope of our model. The 
problem with the jet proposal is that the additional component has to be quite compact 
15 AU) in order not to be seen in the high resolution VLBI images of Marcaide et al. 
(1992) and Alberdi et al. (1993) at 22 GHz. If the low frequency radio emission is from 
the accretion flow itself, then we need to identify a mechanism which could increase the 
radio emission over and above our baseline model. 

In our model, the synchrotron radiation at different frequencies are produced at different 
radii in the flow, the emission at 10^^ Hz coming from close to the black hole and the emission 
at lower frequencies coming from farther out. Thus, one simple way of improving the £t to 
the radio data is to modify the electron temperature outside of a few tens of Schwarzschild 
radii. Figure 5 is only an illustration. The solid line in Fig. 5a corresponds to our baseline 
model, while the dashed line is another model which is identical in all respects except that 
we have arbitrarily set Te = 2 x 10® K over the radius range r = 20 — 1000 (Fig. 5b). This 
ad hoc model fits the data well. 

Is there any reason to think that the electrons might have the profile shown in Fig. 5b? 
The answer is tentatively yes, since there are several radial transport mechanisms which 
could drive the electrons to a nearly isothermal state. First, the long mean free path of the 
electrons can lead to fairly strong radial heat conduction (parallel to held lines). Second, 
synchrotron radiation, which has a thermalizing effect on the electrons (Mahadevan & 
Quataert 1997), can also cause signihcant energy diffusion. Finally, the tangled magnetic 
helds which we assume to be present in the how may move outwards as a result of buoyancy 
and may dissipate their energy and heat the electrons at larger radii. Some of these ehects 
could well be episodic; it is interesting in this connection that the radio hux of Sgr A* is 
known to be quite variable (Zhao et al. 1989). These ehects need further study. 


5. Discussion 

The main result of this paper is that the two-temperature ADAF model (Narayan 
et al. 1995, Rees 1982) provides a viable explanation of the spectral properties of Sgr 



A*. Our basic model, shown by the solid line in Fig. 1, hts most of the high resolution 
measurements from the high frequency (> 86 GHz) radio to the X-ray band, including the 
stringent infrared limit of Menten et al. (1997) and the revised X-ray luminosity. 

The ADAF model described in this paper makes use of a two-temperature plasma, as 
originally described by Shapiro et al. (1976). This involves two key assumptions. First, it 
is assumed that most of the energy released by viscosity goes into the ions and not the elec¬ 
trons. Second, the only coupling between ions and electrons is assumed to be via Coulomb 
collisions. At the present time, these assumptions have not been proved theoretically to 
be valid, nor have they been disproved (e.g. Begelman & Chuieh 1988). Rees et al. (1982) 
made the interesting point that the plasma physics involved is so complex that perhaps the 
best way of testing the two-temperature paradigm is by comparison of astrophysical mod¬ 
els with observations. Fabian & Rees (1995) declared, based on the success of our earlier 
model of Sgr A* (Narayan et al. 1995), that the two-temperature assumption appears to be 
supported by observations. As the present paper shows, the two-temperature ADAF model 
explains the data on Sgr A* even better than it seemed in the 1995 paper. Furthermore, 
the ADAF model has been snccessfnlly applied to several other sonrces (see Narayan 1997 
for a review), making the Fabian & Rees argument even stronger. 

Althongh the ADAF model presented here includes a number of parameters, only one 
parameter is adjusted in the spectral hts, and the qnality of the agreement between the 
model and the data in Fig. 1 must be judged in this context. The free parameter in the 
model is the Eddington-scaled mass accretion rate m which is adjusted so as to ht the 
X-ray luminosity in the ROSAT band. Althongh there is some nncertainty in the X-ray 
Inminosity because of uncertainty in the hydrogen colnmn (cf. § 2 ), we feel this does not 
seriously affect the model ht (see Fig. 4c which shows how small adjustments in the htted 
parameter m will allow the model to remain consistent with any fntnre revisions of the 
X-ray Inminosity). 

Apart from m all the other parameters in the model are assigned standard valnes 
(these valnes are not changed from one application to another). For instance, we invoke 
strict eqnipartition and assign f3 = 0.5. Once f3 is given, the model has a nnique prescription 
(Esin 1997) for the adiabatic index of the gas: 7 = (8 — 3/5)/(6 — 3/3). Farther, by making 
use of the scaling between the shear stress and the magnetic pressure derived by Hawley, 
Gammie & Balbus (1996), we are also able to hx the viscosity parameter: a 0 . 6 ( 1 -/3). 
Of course, the use of the a prescription itself involves an assumption about the nature of 
the viscosity. Narayan (1996), however, has argued that this viscosity prescription with a 
independent of radius is particularly appropriate for ADAFs. The hnal parameter in the 
model is S, the fraction of the viscous energy that goes directly into the electrons. By 
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assumption S is small, and it is usually assigned a value 5 = 10 ^ ~ nie/mp. The results 
are, however, insensitive to the exact value so long as <5 ^ 10“^ (see Fig. 4d). 

The models presented here are substantially more elaborate (cf. §2) than the ones 
described in Narayan et ah (1995). It is pleasing that the more detailed analysis presented 
here leads to an improvement in the agreement with observational constraints. In contrast 
to the models presented in Narayan et al. (1995), where a black hole mass of 7 x IO^Mq 
was required in order to accommodate the radio data, here we hnd that the model quite 
naturally hts the data with a more reasonable mass of 2.5 x IO^Mq (Eckart & Genzel 
1997). This removes one of the main problems with the previous work. The model also 
predicts a mass accretion rate m ~ 1.3 x 10“"^ (in Eddington units) which is within the 
range considered likely by a direct estimate of the accretion rate (Genzel et al. 1994). 

We should emphasize that the model is self-consistent, and includes a detailed treatment 
of hydrodynamics, radiation processes, and thermal balance of ions and electrons in the two- 
temperature plasma. Most other models in the literature are more primitive and/or ad hoc. 
In addition, the model is fairly robust to changes in the parameters (Fig. 4). The predicted 
spectrum is insensitive to large variations in the viscosity parameter a and the electron 
heating parameter S. The results do, however, vary with the equipartition parameter f3 
(Fig. 4b, 4c). 

There are some problems with the present model. First, we under-predict the y-ray flux 
in the EGRET band by an order of magnitude. This might indicate a residual uncertainty 
in the relative normalization of the protons and electrons, but we have not been able 
to come up with any specihc proposal to eliminate the discrepancy. It is worrying that 
we underpredict the flux even though we have assumed a power-law energy distribution 
for the protons which maximizes the predicted y-ray luminosity. Perhaps the y-ray source 
detected by EGRET with its one degree beam is not Sgr A* at all, but diffuse emission 
from a molecular cloud in the Galactic Genter region. 

Another problem is that the model under-predicts the radio flux below 86 GHz. The 
existence of a break in the radio spectrum at around 86 GHz suggests that the region of the 
spectrum below 86 GHz might arise from a distinct component, e.g. a weak outflow of some 
kind, as postulated in the model due to Falcke (1996). We note in this connection that 
Krichbaum et al. (1993) claim to have seen evidence for a jet in their 43 GHz VLBI image, 
though Backer et al. 1993 could not verify this. One problem with the jet hypothesis is 
that the 22 GHz VLBI maps of Marcaide et al. (1992) and Alberdi et al. (1993) show no 
evidence for any resolved structure. This means that the jet has to be smaller than 15 AU, 
i.e. no larger than ten times the size we would predict for the ADAF at this frequency. If 
the radio flux is from a jet, then clearly the discrepancy we face in the radio luminosity is 
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not an issue since our model includes only the accretion flow. However, if the radiation is 
not from a separate component but from the ADAF itself, then it probably implies that 
the electron temperature in the flow differs from that predicted by our model. Figure 5 
shows the kind of temperature profile that is needed to fit the observations. The profile 
shown here is completely ad hoc. However, if radial transport processes are included, the 
model might give a qualitatively similar temperature profile. We would emphasize that the 
luminosity of Sgr A* in low frequency radio waves is a small fraction (less than a percent) 
of the bolometric luminosity of the source. Minor changes in the flow profile can explain 
the discrepancy here without affecting the overall flow dynamics or energetics. 

Yet another problem is the somewhat low electron brightness temperatures predicted 
by the model. As Fig. 4c shows, we obtain values below 10^'^ K at all radii, whereas VLBI 
observations at 43 GHz and 86 GHz indicate brightness temperatures in excess of 10^° 
K (Backer et ah 1993, Rogers et ah 1994). Thus, although the model fits the 86 GHz 
flux well, it seems to predict a larger angular size for the source, by a factor of ~ 1.5, 
than observations indicate. For comparison with future observations, we list here linear 
source radii at various radio frequencies according to the baseline model shown in Fig. 1: 
1.3 X 10^^ cm at 10^^'^ Hz, 4.1 x 10^^ cm at 10^^ Hz, 1.0 x 10^^ cm at 10^^'^ Hz, 2.1 x 10^^ 
cm at 10^^ Hz, 4.1 x 10^^ cm at 10^*^'^ Hz, 8.3 x 10^^ cm at 10^° Hz. Note that the electron 
temperature depends on the value of the parameter (3 (Fig. 4c). Thus, one might be able to 
fit the observed brightness temperatures by tuning this parameter as a function of radius, 
although we have kept (3 fixed in our models. 

An outstanding feature of the ADAF model presented here is its extraordinarily low 
radiative efficiency: e = Lboi/Afc^ = 5 x 10“®. The low efficiency permits the model to fit 
the very low luminosity of Sgr A* with a fairly large m. Figure 1 shows two thin disk 
models with the standard radiative efficiency of e ~ 0.1. A thin disk with the “correct” m 
(dotted line) over-predicts the infrared flux by four or five orders of magnitude. A model 
with m reduced by a factor of 10^ accommodates the IR limits but does not fit any of the 
other data and has a mass accretion rate which is extremely discrepant with independent 
estimates of rh. Some models in the literature which are based on a thin disk attempt 
to solve the luminosity problem by hiding a large fraction of the emission in the optical 
band where there is extinction by dust (e.g. Falcke et ah 1993a, Melia 1994). However, the 
models invariably predict emission in the infrared at a level well above the Menten et al. 
(1997) limit. Other models in the literature which are not based on a thin disk also have 
problems fitting the low luminosity (e.g. Duschl & Lesch 1994; Mastichiadis & Ozernoy 
1994) and generally require values of m significantly below the lower limit of 10“^. 

Gan we save the thin disk model by assuming that the disk has no viscosity at all. 



so that the gas in the disk does not accrete at the present time? In a very interesting 
argument, Falcke & Melia (1997) show that even such an extreme model is inconsistent. 
This is because the mere presence of the disk, even though it does not accrete, will still lead 
to fairly strong infrared emission. The argument is that if there is inflow of material towards 
the Galactic Center via a Bondi-Hoyle-like flow, when the inflowing gas hits the disk and 
circularizes at a radius rdrc it will produce a substantial amount of thermal radiation. For 
the mass accretion rate of m > 10“^ estimated in the case of Sgr A* (cf. §4) the predicted 
infrared flux is well above the Menten et ah (1997) limit in nearly all the cases considered 
by Falcke & Melia (1997). The only way of avoiding this argument is by making rdrc large 
(say Tcirc ^ 10^), but Falcke & Melia show that this is not consistent with most reasonable 
wind specihc angular momentum conhgurations. 

The ADAF model is able to circumvent the Falcke & Melia argument quite naturally 
and this is another point in favor of this model. Since the accreting gas in the ADAF has a 
quasi-spherical shape (Narayan & Yi 1995a), the incoming Bondi-Hoyle flow never reaches 
Tcirc- Instead it shocks with the ADAF at a distance on the order of the outer radius i?out 
(recall that 77 ~ i? in an ADAF). In our models we have taken 7?out = 10^7?*, but as 
we argued in §3 the outer radius may be even larger. The luminosity associated with the 
impact at such large radii is quite low. Once the stream has been incorporated into the 
ADAF at the outer radius, the accretion switches to an advection-dominated form and the 
gas does not radiate very much as it moves in. 

Several recent papers have made the point that the successful application of the ADAF 
model to any observed system is direct evidence that the accreting object is a black hole 
(Narayan, McClintock & Yi 1996; Narayan, Yi & Mahadevan 1996; Narayan, Barret & 
McClintock 1997; Narayan, Garcia & McGlintock 1997). The argument is that if accretion 
is via an ADAF and if the object has an event horizon, then the advected energy will 
disappear from sight. However, if the central object has a surface, then the surface will be 
heated by the hot inflow from the ADAF and the advected thermal energy will be emitted 
as thermal radiation. Thus, for an object with a surface, the radiative efficiency will be 
the canonical 10% even if the accretion occurs via an ADAF. Only if the central object is 
a black hole can the radiative efficiency be truly low. 

Sgr A* is perhaps one of the best objects for this argument. This source provides two 
strong observational constraints: (1) The mass accretion rate is estimated to be at least 
a fewx 10 “®Mq yr“^, and (2) the bolometric luminosity is no greater than 10^'^ ergs“^. In 
our opinion the only plausible resolution of these two conflicting pieces of evidence is to 
postulate (1) that the accretion in Sgr A* occurs via an ADAF and (2) that the central 
object is a supermassive black hole. This “proof” of the black hole nature of Sgr A* is 
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qualitatively different from usual proofs which rely on a measurement of the mass. The 
usual argument is that if an object is too massive to be a neutron star it must be a black 
hole; it is a proof by elimination. Our “proof” is somewhat more direct and cuts to the 
essence of what constitutes a black hole, namely the presence of an event horizon. We argue 
that Sgr A* has an enormous luminosity dehcit for plausible values of M and the simplest 
explanation is that the object has an event horizon. 

One of the results to come out of this work is our conhrmation that compressive heating 
of electrons cannot be ignored and must be included consistently in computations of ADAF 
spectra of low luminosity systems like Sgr A* (see Fig. 2). This point was made recently 
in an important paper by Nakamura et ah (1997). In view of this result it would be useful 
to revisit other low-m systems to which the ADAF model has been applied (e.g. soft X-ray 
transients, Narayan, Barret & McClintock 1997; Hameury et ah 1997; NGC 4258, Lasota et 
ah 1996; low-luminosity nuclei of giant ellipticals, Fabian & Rees 1995; Mahadevan 1997; 
Reynolds et ah 1996) and redo the analysis with the inclusion of compressive heating. 
Interestingly, Esin et al. (1997) have shown that compressive heating has a much weaker 
effect on ADAFs at higher values of m. The reason is that with increasing m the other 
terms in the electron energy equation (1), notably the Coulomb collision term and the 
radiative cooling term become more important, and compressive heating no longer 

dominates. 

In the model presented in this paper the emission in the infrared and the radio (above 
about 86 GHz) comes from fairly small radii ~ lOR*. Since the gas in the ADAF is nearly in 
free-fall, the characteristic time scales of the flow are quite short. We may thus expect rapid 
variability in Sgr A* at these wavelengths. The shortest likely time scale is the dynamical 
time, tdyn (GM/R^) which gives tdyn 5000 s at i? = 10i?s. The viscous time is 
longer than this by a factor of l/a ~ 3. The longer wavelength radio emission is from larger 
radii; for instance, the emission below 1 GHz comes from R ^ lO^i?^. Variability at these 
wavelengths will be correspondingly slower: tdyn(lO^Rs) ~ 5 x 10® s. The bremsstrahlung 
emission in soft X-rays is from a broad range of radii, but is dominated by large radii. 
Therefore, the soft X-ray flux should show much slower variations (timescale ~ 1 year) 
than the radio, millimeter or infrared fluxes. However, the hard X-ray (~ 100 keV) flux at 
the broad peak of the spectrum will trace the higher temperature ( 10^ - lO^® K) 

inner regions and could vary on time scales of days. 

By assumption, the electrons in our model are fully thermal. Mahadevan & Quataert 
(1997) have, however, shown that at low mass accretion rates similar to that in Sgr A* , 
electrons may not be thermalized at larger radii. The electron energy distribution at these 
radii may then be truncated even more sharply than in a Maxwellian, which will cause a 
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reduction in the synchrotron emission. This will act to increase the discrepancy between 
the model and the data at low radio frequencies. This is an area for further work. 

Note, however, that a power-law distribution of electrons extending over any reason¬ 
able range of energies is constrained by the observations. For instance, since optically thin 
nonthermal synchrotron emission usually has a spectral form Ly ~ this means that in 

a vLy plot the optically thin emission would continue to rise as Such a rise beyond 

~ 10^^ — lO^^Hz is ruled out by the observed upper limits in the infrared. 

In this connection, we note that Falcke (1996), Duschl & Lesch (1994) and Beckert & 
Duschl (1997) have proposed nonthermal models for Sgr A* in which most of the electrons 
have Lorentz factors of around a few hundred. These authors require a mono-energetic elec¬ 
tron distribution in order to reproduce the sharp cutoff observed in the sub-millimeter/FIR 
band. Falcke assumes a mono-energetic distribution in the “nozzle” of his jet-disk model, 
while Duschl & Lesch assume a homogeneous sphere of mono-energetic electrons. Such 
mono-energetic models are not ruled out by the argument of the previous paragraph and 
indeed both groups are able to £t the radio spectrum of Sgr A* reasonably well, although 
the sub-mm, IR and X-ray spectrum are not explained. We, however, find the idea of a 
nonthermal but mono-energetic distribution somewhat artificial. Our ADAF model fits the 
data in a more natural way by making use of a thermal distribution of electrons. Fur¬ 
ther, because of thermal synchrotron self-absorption, the cutoff of the spectrum in the 
sub-millimeter band is very sharp without requiring any fine-tuning. Melia’s (1992, 1994) 
model also makes use of self-absorbed thermal synchrotron emission and has similar prop¬ 
erties, but there is an error in his calculation of the synchrotron emissivity (cf. Mahadevan, 
Narayan & Yi 1996). 

Further high spatial resolution observations of Sgr A* are highly desirable. The ADAF 
model can be tested and constrained with better observations of the spectrum and vari¬ 
ability, especially in the millimeter, sub-millimeter, infrared and X-ray bands. Better y-ray 
observations might also resolve the issue of whether the source detected by EGRET is Sgr 
AT 


After this paper was submitted we saw a closely related preprint by Manmoto, Mi- 
neshige & Kusunose (1997) entitled “Spectrum of optically thin advection dominated accre¬ 
tion flow around a black hole: application to Sgr A*.” There are many points of similarity 
between their paper and ours and fairly close agreement in the results. In particular, Man¬ 
moto et ah emphasize the importance of the advection term in the heating of electrons. 
The main difference between the papers is in the choice of equation of state of the electrons. 
We allow for the effect of magnetic fields in calculating the adiabatic index of the electrons 
and we also include relativistic corrections (see Appendix A) whereas Manmoto et al. em- 



ploy a simpler prescription. Also, we make use of a more comprehensive and updated set 
of spectral data for the model comparisons. 


This work was supported in part by NSF grant AST 9423209 and NASA grant NAG 5- 
2837. We thank Shoji Kato for sending us preprints of Nakamura et ah (1997) and Manmoto 
et ah (1997) prior to publication, and T. Krichbaum, J. Marcaide, S. Stolovy and E. Serabyn 
for comments on the spectral data. RM thanks Jun-Hui Zhao and Mark Reid for useful 
discussions on the radio observations of the Galactic Genter. 
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A. Energy Advection by Electrons 


The advection term in equation (1) can be written per unit volume as 

= pTeV^, (Al) 

where Se is the entropy of the electrons per unit total gas mass. This term was ignored 
in much of the earlier work, on the assumption that the temperature of the electrons is 
determined primarily by a balance between the Coulomb transfer term and the radiative 
cooling term in equation (1). However, Nakamura et al. (1997) have shown that energy 
advection by electrons can be important under some circumstances and may even play a 
dominant role in determining the electron temperature. This is the case especially when 
the mass accretion rate is low, as in Sgr A*. 

We begin with the relation 


TedSe = dUe + Pgd 




(A2) 


where Ue is the internal energy of the electrons per unit mass and Pg is the electron pressure. 
We consider a mixture of gas and magnetic helds. If gas pressure contributes a constant 
fraction f3 to the total pressure Ptot, then 


_ pkTi pkTe 
Pot — -1-r 


pkTi ^ pkTe 


IMiTTlu 


Perriu 247r fdpem. 


(A3) 


It seems natural to denote the two terms on the right as the effective ion and electron 
pressure, each including an appropriate fraction of the magnetic pressure. Therefore, we 
have Pe = pkTe/{f3pemu). 

The internal energy of the gas is a sum of the ion, electron and magnetic held internal 
energies: 


3 kTi kTe 

u = - -h a{Te) - 

2 


+ 


6 — 3/5 kTi 
Snp 2(5 PiTTii 


+ 


3(1-/5) 

fd 


+ fl(Te) 




(A4) 


where the coefficient a(Te) varies from 3/2 in the case of a non-relativistic electron gas, to 3 
for fully relativistic electrons. The general expression for a as a function of the dimensionless 
electron temperature 9e = kTe/nieC^ was derived by Chandrasekhar (1939, Chapter X, 
eq.[236]): 

1 /'3A3(1/9.) + A'i(1/9.) 

= T. I— ipm — 


(A5) 
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Note that the ions never become relativistic in these flows, so that the corresponding 
coefficient for the ions is always ~ 3/2. 

As we have done for the pressnre, the right hand side of Eq. (|A4D may be natnrally 
divided into two terms, the ion and electron internal energies: 


6-3/3 kTi 
2[3 fiifriu ’ 


and Ue 


3 ( 1 -/ 3 ) 




kT, 


l^e^u 


(A6) 


In this interpretation, Ui and Ug are again “effective” qnantities, which inclnde contribntions 
from the particles as well as the associated magnetic held. Note that the contribntion of 
the magnetic held to the internal energy of each particle species is proportional to the 
contribntion of these particles to the total pressnre, a natnral choice in onr model where 
the ratio of the magnetic to gas pressnre is fixed. Having defined Pg and Ug, we can now 
write the energy advection term for the electrons as 


Q 


e,adv 


pvk 


3 ( 1 -/ 3 ) 


+ Oj{Tg) + Tg 


da 

dTg 


dR 


vkTg dp 
{dpgrriu dR 


(A7) 


/^From the pressnre and internal energy of the electrons, we calculate the effective 
adiabatic index of the electrons via the relation 7 e — 1 = Pg/ug. This gives 


4 — 3/3 + a/3 
3 — 3/3 + a/3 


(AS) 


If the particles are relativistic, then a = 3 and 7 e = 4/3 regardless of /3. This is because both 
the particles and the tangled held behave like radiation. If the particles are non-relativistic, 
however, a = 3/2 and in this case 


as shown by Esin (1997). We see that 7e varies from 4/3 when /3 = 0 to 5/3 when /3 = 1. 
Thus, the adiabatic index depends fairly sensitively on (3 so long as the electrons are non- 
relativistic. If the electrons in the accreting gas behave adiabatically, then their temperature 
varies as Tg oc (technically there is also the term proportional to da/dTg in eq A7 

which becomes important when the electrons are quasi-relativistic, but we neglect this for 
simplicity). Thus, for a given amount of compression, the electrons are hotter when /3 is 
large and cooler when /3 is small. This explains the trend seen in Figs. 4b and 4c. 
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Figure Captions 

Figure 1. The open and filled circles represent various fiux measurements and upper limits 
of Sgr A*. We consider the filled circles to be more important as model constraints. The 
box at 1 keV represents the uncertainty in the X-ray fiux. The solid line is our baseline 
ADAF model with the following parameters: M = 2.5 x IO^Mq, a = 0.3, f3 = 0.5, 5 = 0.001. 
The Eddington-scaled mass accretion rate m has been adjusted to fit the X-ray fiux, giving 
m = 1.3 X 10“"^. The peak at the left is due to synchrotron radiation, the next peak is 
due to Compton scattering, the peak between 10-100 keV is due to bremsstrahlung and 
the peak above 100 MeV is due to pion production. The long-dashed line is a model in 
which the pion peak has been artificially raised by about an order of magnitude to fit 
the data. The dotted line is the spectrum corresponding to a standard thin accretion disk 
with m = 10“"^ while the short-dashed line is a thin disk with m = 10“^. Neither of these 
models is satisfactory. 

Figure 2. (a) Spectra corresponding to five models, showing the effect of various approx¬ 
imations; see the text for details, (b) Electron temperature profiles corresponding to the 
same five models. The models represented with long-dashed and solid lines have nearly 
identical temperatures. 

Figure 3. (a) The solid line is the baseline model shown in Fig. 1, with a black hole mass of 
M = 2.5 X IO^Mq. The dotted line corresponds to a model with M = 8.3 x IO^Mq and the 
dashed line corresponds to M = 7.5 x IO^Mq. (b) The solid line is the baseline model with a 
mass accretion rate m = 1.3 x 10“"^. The dotted, short-dashed, long-dashed and dash-dotted 
lines corresponds to models with rh = 6.4 x 10“^, 9.0 x 10“^, 1.8 x 10“*^, rh = 2.5 x 10“"^, 
respectively. 

Figure 4. (a) The solid line is the baseline model shown in Fig. 1, with a = 0.3. The 
dotted, short-dashed and long-dashed lines correspond to models with a = 0.1, 0.2, 0.4 
respectively, (b) The solid line is the baseline model shown in Fig. 1, with (3 = 0.5. The 
dotted, short-dashed, long-dashed lines and dot-dashed lines correspond to models with 
f3 = 0.3, 0.4, 0.6, 0.7 respectively, (c) Electron temperature profiles corresponding to the 
models shown in (b). (d) The solid line is the baseline model shown in Fig. 1, with S = 
0.001. The dotted, short-dashed, long-dashed and dot-dashed lines correspond to models 
with (5 = 0, 0.003162, 0.01, 0.03162 respectively. 

Figure 5. (a) The solid line is the baseline model shown in Fig. 1. The dotted line corre¬ 
sponds to a model in which the temperature of the electrons has been fixed at Tg = 2 x 10® 
K over the radius range r = 20 — 1000. This model fits the low frequency radio data better, 
(b) Electron temperature profiles corresponding to the two models shown in (a). 



Table 1: Sgr A*: Radio and NIR Observations (Distance = 8.5 kpc) 


V 

A 

e 

Su 

VLy 

Ref. 

Hz 

/im 

// 

Jy 

ergs s“^ 



Davies et al. 76 


4.08x103 

735294 

9.6x10® 

312500 

9.6x10® 

312500 

9.6x10® 

312500 

1.66x10^ 

180722 

1.5 xlO^ 

200000 

1.5 xlO^ 

200000 

2.3 xlO^ 

130000 

2.7 xlO^ 

110000 

2.7 xlO^ 

110000 

5.0 xlO^ 

60000 

5.0 xlO^ 

60000 

8.1 xlO^ 

37000 

8.1 xlO^ 

37000 

8.3 xlO^ 

36000 

8.4 xlO^ 

35714 

8.4 xlO^ 

35714 

1.5 xlO^o 

20000 

1.5 xlOi° 

20000 


4.3 

< 0.05 

10 

0.29 

10 

0.27 

10 

0.30 

2.5 

0.56 

~ 1.2 

0.8 

~ 1.2 

0.3 

0.2 

1.1±0.1 

~ 0.65 

.73 

~ 0.65 

.42 

~ 0.35 

1.13 

~ 0.35 

0.55 

~ 0.22 

0.9 

~ 0.22 

0.58 

0.016 

1.2±0.1 

~ 0.21 

1.07 

~ 0.21 

0.55 

~ 0.12 

1.64 

~ 0.12 

.68 


< 1.76 X 10^0 
2.41 X 10^^ 
2.24 X lO^i 
2.49 X 10^^ 
8.04 X lO^i 

1.04x10^2 

3.89x10^^ 

(2.2 ±0.2) X 10^2 
1.7 X 10^2 
9.78 X lO^i 

1.45x10^2 

7.13x10^^ 

6.23x10^2 

4.01 X 10^2 
(8.6 ±0.7) X 10^2 
7.77 X 10^2 
3.99 X 1032 
2.1x1033 
8.7x1032 


Davies et al. 76 
Davies et al. 76 
Davies et al. 76 
Davies et al. 76 
Backer 82 
Zhao et al. 89 
Marcaide et al. 92 
Backer 82 
Brown & Lo 82 
Zhao et al. 89 
Zhao et al. 89 
Backer 82 
Brown & Lo 82 
Marcaide et al. 92 
Zhao et al. 92 
Zhao et al. 92 
Zhao et al. 92 
Zhao et al. 92 






Table 1: Sgr A*: Radio and NIR Observations (continued). 


V 

A 

d 

Su 

VLy 

Ref. 

Hz 

/xm 

// 

Jy 

ergs s“^ 



1.5 xlO^o 20000 

2.2 xl0^° 13600 

2.2 xl0^° 13600 

2.2x10^° 13500 

2.3x10^° 13000 

x4.3xl0i° 7000 

8.6 xlO^o 3488 

8.6 xlO^o 3488 

X 8.6 xlO^o 3488 
2.2x10^^ 1350 

2.3x10^^ 1300 

2.3x10^^ 1300 

3.5x10^1 870 

3.75x10^^ 800 

5.0x10^1 600 

6.7x10^^ 450 

6.7x10^1 450 

8.6x10^1 350 

8.6x10^1 350 

8.6x10^1 350 

1.0x10^3 30 


0.12x0.24 

1.15+0.01 

~ 0.08 

2.1 

~ 0.08 

0.8 

1.8x10-3 

1.07 + 0.15 

1.8x10-3 

1.2 ±0.4 

0.75x10-3 

1.4+0.1 

0.02 

1.3 

4x8 

1.05 

0.16 X 10-3 

1.4+0.2 

1.9x4.3 

2.4 ±0.5 

11 

2.5 

11 

2.6 ±0.6 

8 

4.8 + 1.2 

13 

3.5 ±0.5 

10 

4.0 + 1.2 

7 

< 1.5 

8 

3.0 + 1.0 

11 

7 + 2 

11 

< 10 

30 

18.5 + 9 

8 

< 120 


(1.49+0.06) X 10^3 

3.99 X 1033 
1.52x1033 

(2.04 + 0.29) X 1033 
(2.4 ±0.8) X 1033 
(5.2 ± 0.4) X 1033 
9.67x1033 
7.8x1033 

(1.04 + 0.15) X 1034 
(4.57+ .95) X 1034 
4.97 X 1034 
(5.2 + 1.2) X 1034 
(1.45 + 0.36) X 1033 
(1.13 + 0.16) X 1033 
(5.19 ± 0.52) X 1033 

< 8.69 X 1034 
(1.74 + 0.58) X 1033 

(5.2 + 1.5) X 1033 

< 7.4 X 1033 
(1.38 + 0.67) X 1036 

< 1.04 X 1033 


Yusef-Zadeh et al. 90 
Zhao et al. 92 
Zhao et al. 92 
Alberdi et al. 93 
Marcaide et al. 92 
Krichbaum et al. 94 
Backer 82 
Wright et al. 87 
Rogers et al. 94 
Serabyn et al. 92 
Zylka & Mezger 88 
Zylka et al. 92 
Zylka et al. 92 
Zylka et al. 95 
Zylka et al. 95 
Dent et al. 93 
Zylka et al. 95 
Serabyn et al. 97 
Mezger 1994 
Zylka et al. 92 
Zylka et al. 92 



Table 1: Sgr A*: Radio and NIR Observations (continued). 


V 

A 

d 

s. 


Ref. 

Hz 

/iUl 

// 

Jy 

ergs s”! 


1.0x10^3 

30 

4 

< 20 

< 1.7 X 1037 

Telesco et al. 96 

X 1.5x10^3 

20 

1.6 

< 1 

< 1.3 X 1036 

Gezari et al. 94 

1.56x1013 

19.2 

4 

< 1.4 

< 1.9 X 1036 

Telesco et al. 96 

X 1.7x1013 

18 

2.3 xl.3 

< 0.3 

< 4.4 X 1036 

in Zylka et al. 92 

2.3-3.6x1013 

13-8 

2.3 xl.3 

< 0.1 

< 2.6 X 1036 

in Zylka et al. 92 

X 3.5x1013 

8.7 

0.5 

< 0.1 

< 3.0 X 1036 

Stolovy et al. 96 

X 1.4x101^ 

2.2 

0.15 

< 9 X 10-3 

< 1.1 X 1036 

Menten et al. 97 



Table 2: Sgr A*: X-Ray and 7 -Ray Observations (Distance = 8.5 kpc). 


Energy 

Telescope/ 

e 

Leb ^ 

VLy 

Ref. 

Band (EB) 

Instrument 


erg s“^ 

erg s“3 


0.8 - 2.5 keV 

ROSAT 

~ 20 " 

1.55 X 1034 

1.6 xl034 

Predehl & Trumper 94 ^ 

2 - 10 keV 

ASGA 


< 6.4 X 1033 

< 4.8 X 1035 

Koyama et al. 96 

35 - 75 keV 

SIGMA 

~ 15' 

< 3.5 X 1035 

< 4.8 X 1035 

Goldwrum et al. 94 

75 - 150 keV 

SIGMA 

~ 15' 

< 2.4 X 1035 

< 3.6 X 1035 

Goldwrum et al. 94 

30 - 50 MeV 

EGRET 

~ 1° 


< 1.8 X 1035 

Merck et al. 96 

50 - 70 MeV 

EGRET 

1 0 
~ 1 


< 1.1 X 1035 

Merck et al. 96 

70 - 100 MeV 

EGRET 

1 0 
~ 1 


< 8.3 X 1035 

Merck et al. 96 

100 - 150 MeV 

EGRET 

1 0 
~ 1 


< 6.2 X 1035 

Merck et al. 96 

150 - 300 MeV 

EGRET 

1 0 
~ 1 


(4.9li-?) X 1035 

Merck et al. 96 

300 - 500 MeV 

EGRET 

1 0 
~ 1 


(1.2l°J^) X 1035 

Merck et al. 96 

500 - 1000 MeV 

EGRET 

1 0 
~ 1 


(1.415:1) X 1035 

Merck et al. 96 

1 - 2 GeV 

EGRET 

1 0 
~ 1 


(1.7+5f) X 1035 

Merck et al. 96 

2 - 4 GeV 

EGRET 

1 0 
~ 1 


(2.2l5:i) X 1035 

Merck et al. 96 

4 - 10 GeV 

EGRET 

1 0 
~ 1 


(8.3lt?) X 1035 

Merck et al. 96 


is the total luminosity integrated over the band. 

^ This flux is obtained using Nh = 6 x as opposed to the much higher 

column density used by Predehl & Trumper (1994). This is discussed in §2. 









log [T^(K)] log[yL^(erg s-‘)] 



















